Comment on "Critical behavior of the Pauli spin susceptibility, 
by A. A. Shashkin et al., |cond-mat/0409100 
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The paper by A. A. Shashkin et al. reports measurements of the thermodynamic magnetization of 
two-dimensional electrons in silicon. Although the experimental data is very similar to that reported 
by us 0] more then two years ago, the authors arrive at an opposite conclusion regarding the spin 
susceptibility "critical behavior" and spin instability in the vicinity of the metal-insulator transition. 
We show that this interpretation is based on a flawed analysis of the experimental data. 



The state of strongly correlated fermions at low tem- 
peratures, and in particular its magnetic properties, that 
has been a subject of active interest, is still poorly un- 
derstood. Some theories predict that in the absence of 
disorder exchange effects can lead to spontaneous mag- 
netization, but such instability has never been observed, 
see for references. Several years ago we proposed a 
novel technique for measuring the thermodynamic mag- 
netization of 2d electrons and applied it to a high mo- 
bility silicon MOSFET. Notwithstanding the enhanced 
spin susceptibility observed in the experiment, the mag- 
netization vanished smoothly as the magnetic field was 
reduced, thus excluding the possibility of a spin instabil- 
ity at zero field. 

Ref. 0] repeats our measurements and reproduces the 
data, extending it to 20% lower carrier concentration. Al- 
beit the similarity in the data, the conclusion of Ref. Q 
is opposite to ours. Based on extrapolation from strong 
fields authors infer magnetic instability at zero field. 
In what follows we show how this erroneous conclusion 
emerges from a wrong assumption. 

The experimental method is described in Prus et al. . 
It yields directly dfi/dB which by the Maxwell relation 
equals —dM/dn, with /i, n the electron chemical poten- 
tial and density, and M, B the magnetization and mag- 
netic field. Fig. 2 in Ref. Q presents d^/dB for B be- 
tween 1.5 T and 7 T. The curves @ for B = 1.5, 7T are 
reproduced as dM/dn by solid lines in Fig.^ In Ref. |2j 
the vanishing of du/dB is identified with the full spin 
polarization. The corresponding B values plotted versus 
density (Fig. 4a in Ref. 0) align along a straight line 
that, when extrapolated to B — 0, crosses the density 
axis at a finite n close to the critical density n c (0) below 
which the system becomes strongly localized at B = 0. 
From such an extrapolation Ref. |2j concludes that the 
crossing indicates spontaneous spin polarization at zero 
field. 

We first note that the association of du/dB = 
with full spin polarization is baseless. According to the 
Maxwell relation, full spin polarization would occur at 
du/dB = —dM/dn = —ub, where /j,b is the Bohr mag- 
neton. In contrast, dM/dn — occurs at the maximal 
magnetization rather then at the full one, and we denote 
the corresponding density n m . Fig. Et: taken from 0, 
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FIG. 1: Solid line: dM/dn replotted from 0. Dashed line: 
the same data shifted by 0.15/is and 0.7/is for B — 1.5 and 
7T respectively in order to account for the diamagnetic con- 
tribution of — 0.1/xb/ T. 



depicts spin magnetization vs. density obtained by inte- 
gration of dM/dn data. We note that the density n m (B) 
is always smaller than n c (B), and the magnetization at 
this density, M(B,n m (B)), is always below the full one, 
HBn m (B). 

The degree to which the assumption of full spin po- 
larization at n m is invalid is also evident directly from 
the data of Ref. ;2]. The dependence M(n rn ) can be ob- 

tained by integration: M(n m ) — J dndM/dn. Since 

o 

the measured dM/dn (Fig.^l never approaches /j,b, the 
maximal magnetization is significantly lower then the full 
one, and therefore the susceptibility determined in [2| as 
HBn m {B)/ B and plotted in fig. 4b of [2J is also signifi- 
cantly overestimated. 

We also remark on the risk of drawing conclusions at 
-6 = based on extrapolation from high B. It is quite 
easy to fall into the trap of arguing for instability at low 
B by attempting to judge on the faith of the maximal 
magnetization or the magnetization at n c by the behavior 
at high field. However, as clearly seen in Fig. 03, both 
go smoothly to zero with vanishing magnetic field. 
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FIG. 2: (a) Spin magnetization as a function of density at 
different magnetic fields and temperatures 0.2, 0.8, 2.5 and 
4.2 K; higher magnetization corresponds to lower tempera- 
ture. Critical densities, n c , are marked by circles. Thick blue 
line - full magnetization, thick red line - magnetization of a 
degenerate ideal electron gas at B = 6 T. (b) Maximal spin 
magnetization and spin magnetization at the critical densities 
plotted against magnetic field. Dashed line - extrapolation 
from high magnetic fields (reproduced from Q). 



Finally we note that one should not ignore the diamag- 



netic contribution to M 3]. This contribution can be 
roughly estimated from the experimental data by assum- 
ing that the maximal magnetization at high fields, e.g. 
7 T, approaches lfis per electron. To obtain such value 
from Fig.^the curve should be shifted by approximately 
0.7/jIb- This will also increase dM/dn for 1.5 T, keep- 
ing it still well below the full spin polarization value (see 
dashed lines in Fig. This approximation for the dia- 
magnetic shift is supported by theoretical calculations Q] . 

In conclusion, the data presented in Ref. @|, properly 
analyzed, indicates that the 20% reduction in the accessi- 
ble density compared to the earlier studied samples does 
not lead to spin instability. Such an instability, antici- 
pated by theory, may still occur at a much lower density. 
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[3] For parabolic confining potential U{x) — m' 'lo 2 x 2 /2, in 
the presence of parallel magnetic field with the cyclotron 
frequency lo c , the electron energy is Eo = h%(u)Q + lj 2 ) 1 ^ 2 . 
Expanding (uq+^) 1/2 = lu + lu 2 / (2lu ) , with duj c /dB = 
e/m*c, the ground state energy field dependence is 

dEo/dB « (m/m*) x (w c /wo) X /is 

The diamagnetic contribution to the energy of spatial 
quantization is linear in B, as expected. For the chemi- 
cal potential set by Eo we obtain d^i/dB w dEo/dB. The 
ratio ui c /ujo can be estimated as (d/X) 2 , with d the 2DEG 
thickness and A the magnetic length. For d w 30A and 
B = 1 T the diamagnetic contribution is comparable to 
\ib and thus cannot be ignored. 



